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Families of stable and metastable solitons in coupled system of scalar fields 
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In this paper, we obtain stable and metastable soliton solutions of a coupled system of 
two real scalar fields with five five discrete points of vacua. These solutions have definite 
topological charges and rest energies and show classical dynamical stability. From a 
quantum point of view, however, the V-type solutions are expected to be unstable and 
decay to D-type solutions. The induced decay of a 1^-type soliton into two D-type ones is 



calculated numerically, and shown to be chiral, in the sense that the decay products do not 

_C ■ respect left-right symmetry 
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INTRODUCTION 



Relativistic solitons, including those of the well-known Sine-Gordon (SG) equation, exhibit re- 
markable similarities with classical particles. They exert short range forces on each other and 
collide, without losing their identities lj-|a|. They are localized and do not disperse while prop- 
agating in the medium. Because of their field nature, they do tunnel a barrier in certain cases, 
although this tunnelling is different from the well-known quantum version [2|, |6[ .Topological soli- 
tons are stable, due to the boundary conditions at spatial infinity. Their existence, therefore, is 
essentially dependent on the presence of degenerate vacua [2, lZ|- 

Topology provides an elegant way of classifying solitons in various sectors according to the 
mappings between the degenerate vacua of the field and the points at spatial infinity. For the Sine- 
Gordon system in 1 + 1 dimensions, these mappings are between (j) = 2n7r, n £ Z and x = ±oo, 
which correspond to kinks and antikinks of the SG system. More complicated mappings occur 
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a. 



in solitons in higher dimensions [3|. Coupled systems of scalar fields with soliton solutions have 



found interesting applications in double-strand, long molecules like the DNA molecules 
bi-dimensional QCD 121 ]. 



HI] and 



Bazeia et al. |13l ] considered a system of two coupled real scalar fields with a particular self- 
interaction potential such that the static solutions are derivable from first order coupled differential 
equations. Riazi et al. [14J employed the same method to investigate the stability of the single- 
soliton solutions of a particular system of this type. Inspired by the coupled system introduced 
in [2| , we propose a new coupled system of two real scalar fields which shows interesting types of 
solitons with well-defined topological charges and rest energies (masses). 

In this paper, we focus on a system with the self-interaction potential 

v(cf>, VO = ^ 2 (V 2 - V> 2 ) 2 + V> V - <Ao) 2 , (i) 

in which eft and ip are real scalar fields, and (fio and ipo are constants. This potential is plotted in 
Fig. ?? for (f)Q = 1 and ipQ = 2. Following the common terminology in field theory, this potential 
has absolute degenerate minima at <j) = tp = 0, <j) = ±</>o and ip = ±ipo known as the true vacua. 
In our proposed system, a spectrum of solitons with different rest energies exists which are stable 
or meta-stable, depending on their energies and boundary conditions. 

The structure of this paper is as follows: in Section [II] we review some basic properties of the 
proposed system. In Section IIIIl some exact solutions together with the corresponding charges 
and energies are derived. The necessary nomenclature and general behavior of the solutions due 
to the boundary conditions are also introduced in this section. Numerical solutions corresponding 
to different boundary conditions are presented, and properties of these solutions like their charges, 
masses, and stability status are addressed in this section. In order to investigate the stability of the 
numerical solutions, their evolution is worked out numerically. The classical dynamical stability of 
the solutions is investigated further in llVl Our conclusions and a summary of the results are given 
in the last section |V] 

II. DYNAMICAL EQUATIONS, CONSERVED CURRENTS, AND TOPOLOGICAL 

CHARGES 

Within a relativistic formulation, the Lagrangian density of the system is given by: 

£ = \d^d^ + \d^d^ - {<f> 2 ^ 2 ~ V> 2 ) 2 + W - <Pl) 2 ]. (2) 





FIG. 1: Left: The self-interaction potential (1) is shown as a height diagram over the (<j),ip) plane. (4>o,i/jq) — 
(1,2) is assumed throughout this paper. Right: The self-interaction potential as a contour map over the 
(0, -0) plane. 

From this Lagrangian density, we obtain the following equations for <j) and ip, respectively: 



□0 = -2«K^-V>oT-W(<r-</>o) 



(3) 



and 



□V> = -2VW - 4>ZY - #<w - r ) 



(4) 



Since the lag rangian density Eq.([2|) is Lorentz invariant, the corresponding energy-momentum 



tensor 



ia 



ig is: 



Tfj, u = d^dycf) + d^dyi) - g^ u C; 



(5) 



which satisfies the conservation law 



dfjr^ = o. 



(6) 



In Equation ©, g^ = diag(l, —1) is the metric of the (1 + 1) dimensional Minkowski spacetime. 
The Hamiltonian (energy) density is obtained from Eq.([5|) according to 



U = T 



oo 



US) 2 +Uf) 2 +HiJ +ns) +^>- 



1 1 fd^ 2 



(7) 



It can be shown that the following topological currents can be defined, which are conserved inde- 
pendently, and lead to quantized charges: 






1 
2^o' 



4 = —^d v ^. 



(8) 



The subscripts H, V and D (to be used later), denote "horizontal", "vertical" and "diagonal" which 
will be explained later. The currents Jjjy are conserved, independent of each other: 



d»r H = 0, 
s M 4 = 0. 



0) 



It is obvious that the total horizontal and vertical charges are conserved separately in any local 
dynamical evolution of the system. The corresponding topological charges are given by: 



/+00 1 

J^cte = -[cK+oo)-0(-oo)], 
-00 ^ 

/+00 -^ 

J^dx = -[^(+00) - V(-oo)]. 



(10) 










FIG. 2: Nomenclature of horizontal (H), vertical (V) and diagonal (D) solutions according to the boundary 
conditions on the ((/>, tp) plane. 

III. SINGLE-SOLITON SOLUTIONS 



Static solutions which correspond to transitions between neighboring vacua are symbolically 
shown in Fig. [2j Accordingly, we call the static solutions H (horizontal), V (vertical) and D 



(diagonal) types, depending on the beginning and end points of the solution on the (<f>, ip) plane. 
We have called them "horizontal" , "vertical" and "Diagonal" simply because of their orientation 
in the (</> ,ip) plane. If we apply -in turn- the constraints ip = tpQ, <p = <Poi an< i "P = Ir*^ i n the 
Lagrangian density ([2]) , the resulting one degree of freedom equations will reduce to the well-known 
<p and p e systems which possess the following exact single-soliton solutions 17], [l8| , summarized 



in Table 1. Of course, these are not the exact solutions of the coupled equations © and 



Table 1. Exact static solutions and the corresponding horizontal and vertical charges. 



Type 



H 

V 
D 



Solution 



<p = ± tanh(2 V2x) , ip Q = ±2 
00 = ±1, ip = ±2 tanh(2\/2x) 

= ±(1/2)V>, V = ±- 



'(l+cxp(±4y / 2g)) 1 / 2 



Qh 



±1 
o 

xl 



Qv 





±1 



Since the functions given in Table 1 are not exact solutions of the coupled equations ([3]) and @ , 
we use them as initial guesses to find the minimum-energy solutions. Any static solution has an 
orbit f(<p, ip) = in the (0, ip) plane, which begins and ends at one of the true vacua. We therefore 
have 



d[_ _ df_ , df_ , 

dx d<p dip 



(11) 



where prime means differentiation with respect to x. The field equations ([3E]) can be integrated 
once to yield 



1 



'\2 



dV 



+ C U 



(12) 



and 



^ f = f^ + c 2 , 



(13) 



where C\ and C2 are integration constants. Combining equations (|lip -(|13 p . we obtain the following 
integro-differential equation for the orbit: 



.Q 



dV 



+ C X 



9/ 

dtp 



dV 

dtp 



dip + Ct 



(14) 



Bazeia et al.[19|| have shown that if the potential can be written in the form 
TT 1/ JF dG\ 2 1/ BF 8G^ 2 



(15) 



then any solution to the following first order coupled equations 



# ^ BF BG n 



and 



<ty dG ,5F =Q (m 

dx dip dtp 

solves the corresponding second order equations. For the coupled system considered in this paper, 
we could use this method to obtain the diagonal (D-type) solutions, provided that 4>q = ipQ. In 
this case, the F and G functions will be 

F = -\\a 2 <P, (18) 

and 

G = iA(0 2 -a 2 )V, (19) 

where a = 4>o- For the D-type solutions, we have (p 2 = tp 2 , and the exact solution will be 

(j) 2 = if: 2 = - a 2 [l + tanh(Aa 2 (x - x ))] . (20) 

We could use this method to find only the D-type solutions. In order to find other stable static 
solutions, we have used the expressions depicted in Table 1 as the starting "guesses" in a step- wise 
variational procedure [2(|. In this procedure, the bi-dimensional spacetime is represented by a grid of 
spatial and temporal step-sizes e and delta, respectively. All the derivatives in the energy density 
expression are then written in their discrete form (e.g. dcp/dx — > ((pi — 0^_i)/e). The initial guess is 
then repeatedly varied in steps, with the total energy checked at each step. Only variations which 
reduce the total energy are accepted and the rest are rejected. The fields are thus varied toward 
the minimum energy solution, were small variations no longer lead to a decrease in total energy. 
The minimum energy, static solutions obtained in this way are plotted in Figure 3. In Figure 4, 
we have mapped these solutions into the (4>, ip) plane. 

The conservation of V and H type charges (Equation ED does not come from a continuous 
symmetry of the Lagrangian, but rather it stems from the boundary conditions and topological 
considerations [20(]. The system under consideration respects various types of symmetries: 1) 
Lorentz symmetry: this symmetry comes from the invariance of the Lagrangian ([2]) under Lorentz 
transformations. Since cp and ip are scalar fields, we conclude that if ((p(x), ip{x)) is a static solution, 
then {cp{^{x — vt),ip( , y(x — vt)) where v is a constant (soliton velocity) and 7 = (1 — v 2 )~ 1 ' 2 , is also 



solution. 2) If (4>(x, t),ip(x, £)) is a solution, then (±0(x, £), ±^(x, £)) are solutions, also. 3) For the 
special case </>o = ipo, we have dual solutions (/> o ip. 4) Parity and time reversal are symmetries of 
the system. Therefore if (4>(x,t),ip(x,t)) is a solution, then (cfi(±x, ±t),ip(±x, ±t)) are solutions, 
too. 
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FIG. 3: Minimu m energy, static solutions. (a) Hbc, (b) Hde, (c) Vd_b, (d) Vec, (e) -Dac and (f) Dea- 
The solid curves represent ip and the dash-dotted curves are for ip. 



IV. SOLITON STABILITY AND CHIRAL DECAY OF F-TYPE SOLITONS 



Basic properties of the H, V, and D solitons are summarized in Table 1. It can be seen that 
all D solitons are degenerate (i.e. have the same rest energy). This degeneracy results from the 
(4>,ip) O (±4>,±ip) symmetry of the Lagrangian density (2). 

Another interesting property of these solitons is that the mass (rest energy) of the Vec soliton 
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FIG. 4: The minimum energy, static solutions mapped on the (</>,V>) plane: (a) Hue, Vdb, Hbc and Vec 
and (b) D DA , D AC , D EA and D AB . 



is larger than the sum of the Dea and Dac solitons, while the topological charge of Vec is equal 
to the sum of the Dea and Dac charges. The same is for Vdb which is more massive than 
the sum of Dda plus Dab- The decay of Vec into Dea + Dac (or Vdb into Dda + Dab) is 
therefore possible with regard to energy and charge conservations. These decays were not observed 
to occur spontaneously in our dynamical simulations, implying that Vec an d Vdb are classically 
metastable configurations. Quantum tunnelling or external perturbations can -in principle- cause 
such a decay (see 2]). In order to show that stimulated decay is in fact possible in the system 
under consideration, we have performed the following simulation: The static vertical soliton Vdb is 
used as the initial condition in a program which calculates the dynamics of the system, by solving 
the coupled PDEs [3] and SJ This soliton is excited, by pumping some energy into the <f> field, via 
an increase in its amplitude. The system is then observed to become unstable and decay into Dda 
and Dab ■ This simulation is shown in Figure Note that the excess energy is transferred into 
both the kinetic energy of the decay products and also emission of some low amplitude waves. 

It is interesting to note that the decay of y-type solitons is chiral, in the sense that when Vec 
decays into Dea and Dac, the former always moves to the left, while the latter always move to the 
right. The fact that the reverse never happens is due to the boundary conditions and topological 
reasons. Similarly, for the decay Vdb — > Dda + Dab, the decay product Dda always moves to 
the left and Dab always moves to the right. This property can be extended to an arbitrary array 
of solitons in a multi-soliton system. The consistency of boundary conditions between successive 
solitons permits certain sequences and forbids certain ones. Examples of allowed and forbidden 



arrays are the followings: 



Allowed arrays : ...D da D ab H B cVceHed- 



Forbidden arrays : ...D da D ba H bc VecHec- 
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FIG. 5: Stimulated decay of a Vdb soliton into Doa + -Dab- Note that the excess energy is transferred into 
the kinetic energy of daughter solitons and also emission of low amplitude waves. Also note that the decay 
is chiral, in the sense that Dda always moves to the left. 
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Table 2. General properties of different soliton families. Note that the soliton masses and 
stability status depend on the choice of (fro and ipQ. Here, we have assumed <f)Q = 1 and tpQ = 2. 



Symbol 


Mass 


Qh 


Qv 


Stability 


Decay Mode 


Dab 


2.858 


-1/2 


+1/2 


stable 


- 


Dac 


2.858 


+1/2 


+1/2 


stable 


- 


D da 


2.858 


+1/2 


+1/2 


stable 


- 


Dea 


2.858 


-1/2 


+1/2 


stable 


- 


Hbc 


3.594 


+1 





stable 


- 


Vec 


6.383 





+1 


metastable 


Vec -> d ea + D A c 


Hde 


3.594 


+1 





stable 


- 


V db 


6.383 





+1 


metastable 


Vdb — > D D a + Dab 



V. CONCLUDING REMARKS 

We studied a nonlinear system of coupled scalar fields which had three (H, V, and D) types 
of stable and metastable soliton solutions. We started by deriving analytical, static solutions 
with definite topological charges. These static solutions, however, were shown to be unstable 
and the stable solutions were obtained numerically, by minimizing the total energy, using a step- 
wise variational method. The full, dynamical equations were used to ensure the stability of these 
minimum energy solutions. According to energy and topological charge conservations, the V-type 
solutions are expected to decay into D-type ones. This, however, does not happen spontaneously 
at the classical level. Either quantum tunnelling or external triggers are needed for decays like 
Vec ~^ Dea + DAC. Conversely, the formation of Vec from Dea + Dac collision is endothermic, 
requiring -at least- the incoming kinetic energy AE = 6.38 — 2 x 2.86 = 0.66 in dimensionless units. 
The stimulated decay of a U-type soliton were calculated numerically, by pumping some energy 
into the (/>-field. This decay is exothermic and the excess energy was observed to be transferred 
into the kinetic energy of decay products and emission of low amplitude waves. 
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